arXiv:1508.05240vl [astro-ph.SR] 21 Aug 2015 


Mon. Not. R. Astron. Soc. 000. [T mH f9999') Printed 24 August 2015 (MN style file v2.2) 


ALMA reveals sunburn: CO dissociation around AGB 
stars in the globular cluster 47 Tucanae 

I. McDonald^*, A. A. Zijlstra^, E. Lagadec^’^, G. C. Sloan^, M. L. Boyer^, 

M. Matsuura^’®, R. J. Smith^, C. L. Smith^, J. A. Yates^, J. Th. van Loon^, 

O. C. Jones^’^, S. Ramstedt^, A. Avison^, K. Justtanont^®, H. Olofsson^*^, 

J. A. D. L. Blommaert^\ S. R. Goldman^, M. A. T. Groenewegen^^ 

^ Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, The University of Manchester, 
Oxford Road, Manchester, M13 9PL, UK 

^Observational Cosmology Lab, Code 665, NASA Goddard Flight Center, Greenbelt, MD 20771, USA 
® The Cornell Center for Astrophysics and Planetary Science, Cornell University, Ithaca, NY If853-6801, USA 
^ Observatoire de la Cote d’Azur, Boulevard de I’Observatoire, CS 34229, F 06304 Nice Cedex 4, France 
^Department of Physics and Astronomy, University College London, Gower Street, London, WCIE 6BT, UK 
^School of Physics and Astronomy, Cardiff University, Queen’s Buildings, The Parade, Cardiff, CF24 3AA, UK 
’^Astrophysics Group, Lennard-Jones Laboratories, Keele University, STS 5BG, UK 
^STScI, 3700 San Martin Drive, Baltimore, MD 21218, USA 

^Department of Physics and Astronomy, Uppsala University, Box 516, SE-75120 Uppsala, Sweden 

Chalmers University of Technology, Dept, of Earth & Space Science, Onsala Space Observatory, S-439 92 Onsala, Sweden 
Astronomy and Astrophysics Research Group, Department of Physics and Astrophysics, Vrije Universiteit Brussel, Pleinlaan 2, 
1050 Brussels, Belgium 

Koninklijke Sterrenwacht van Belgie, Ringlaan 3, B-1180 Brussels, Belgium 


Accepted 9999 December 32. Received 9999 December 32; in original form 9999 December 32 


ABSTRACT 

ALMA observations show a non-detection of carbon monoxide around the four 
most luminous asymptotic giant branch (AGB) stars in the globular cluster 47 Tu¬ 
canae. Stellar evolution models and star counts show that the mass-loss rates from 
these stars should be ~1.2-3.5 x 10“^ Mq yr“^. We would naively expect such stars 
to be detectable at this distance (4.5 kpc). By modelling the ultraviolet radiation 
field from post-AGB stars and white dwarfs in 47 Tuc, we conclude CO should be 
dissociated abnormally close to the stars. We estimate that the CO envelopes will 
be truncated at a few hundred stellar radii from their host stars and that the line 
intensities are about two orders of magnitude below our current detection limits. The 
truncation of CO envelopes should be important for AGB stars in dense clusters. Ob¬ 
serving the CO (3"2) and higher transitions and targeting stars far from the centres of 
clusters should result in the detections needed to measure the outflow velocities from 
these stars. 

Key words: stars: mass-loss — circumstellar matter — infrared: stars — stars: winds, 
outflows — globular clusters: individual: NGC 104 — stars: AGB and post-AGB 


1 INTRODUCTION 

Our understanding of how baryonic matter is recycled from 
metal-poor stars back into their host environment currently 
faces two major problems: how is mass actually lost from 
the stars and what happens to it when it is returned to the 
interstellar medium (ISM)? 


* E-mail: mcdonald@jb.man.ac.uk 


Mass loss from giant stars follows a two-stage pro¬ 
cess. Chromospheric or magneto-acoustically-driven mass 
loss occurs during a star’s red giant branch (RGB) and 
early asymptotic giant branch (AGB) evolution, when 
warm (^^6000 K) plasma is ejected from the surface (e.g . 


Dugree^^^^ 198j; Lim et aP 1 19981 : iLobel fc Dupred 120001 : 


Schroder fc Cunt3l2005l l. Although initially faster, this wind 


slows to ~10 km s at th e luminosity of the RGB tip 
llMcDonald fc van Loon|[2007l : [&oeneweg^l2014l l. Later, n- 
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Table 1. Published properties of the observed stars. 


ID 

RA 

(J2000) 

Dec 

(J2000) 

r 

(') 

< ^LSR > 
(km s“^) 

5v 

(km s“^) 

<5y 

(mag) 

P 

(days) 

L 

(Lo) 

Teff 

(K) 


Notes 

(1) 

(1) 

(2) 

(3) 

(3) 

(3) 

(3) 

(4) (5) 

(4) (5) 

(5) 

VI 

00 24 12.65 - 

-72 06 39.9 

1.87 

-14.7 

10 

4.03 

221 

4824 4760 

3623 3410 

0.90 

V2 

00 24 18.56 - 

-72 07 59.0 

3.26 

-12.7 

12 

2.78 

203 

3031 4470 

3738 3620 

0.35 

V3 

00 25 15.96 - 

-72 03 54.8 

5.49 

-35.2 

11.5 

4.15 

192 

2975 4590 

3153 3540 

0.18 

V8 

00 24 08.59 - 

-72 03 54.9 

0.99 

-32.7 

8 

1.6 

155 

3583 

3578 


References: (l)|Cutri et al.l (1200,'^. 2MASS); (2) 2MASS offset from cluster core 

(00 24 05.67,-72 04 52.6); 

(3) average and 

semi-amplitude of radial velocities fromiLebzelter & WoodI ll2005^: iLebzelter et abl ll2005f) with a —8.7 km s ^ correction 

from heliocentric velocity, plus V-band magnitude pulsation amplitudes Rull range) and periods: f4) [McDonald et al.[ 

(|2011a|). based on multi- 

-wavelength literature photometry: fS) [Lebzelter et ah! ll20l4). based on instantaneous measure- 


ments at the given pulsation phase. 










Table 2. Predicted properties of the observed stellar winds. 


ID 

Mir 

-.dust 

^exp 

'pm.ax 

Dust 

AfReimers 

Msecs 


(Mo yr-i) 

km s~^ 

(K) 

composition 

(Mq yr’ 


(Mo yr-i) 

Notes 

(1) 

(1,3) 

(1) 

(1) 

(2) 


(2) 

VI 

2.1 X 10-® 

4.0 

900 

Silicate, metallic iron 

2-9tt.6 X 

10-1 

Lit®;® X 10-® 

V2 

1.2 X 10“® 

3.8 

900 

Silicate, metallic iron 

2.ltl l 

10-1 

6^4 X 10-1 

V3 

0.9 X 10“® 

3.2 

1000 

Metallic iron only 

2At\l X 

10-1 

7^4 X 10-1 

V8 

1.5 X 10“® 

4.0 

900 

Silicate, alumina, metallic iron, oxides 

-*-•^-0.5 ^ 

10-1 

l.oto 3 ^ 39"® 


References: (1) mass-loss rates, wind expansion vel ocities, and dust condens ation temper atures and mineralogies, base d 
on modelling dust emission in infrared spectr a. from | M cDonal d et al.l ll2011alb (2) based on [McDonald &; Ziilstral ll2015blb 
minimum expected mass-loss rates following iReimer j lll975h with rt = 0.477± 0.070'*~Q'QgS and lSchroder fc Cunt j ||2005| . 
hereafter SC05) with rj = 0.172 ± 0.024^p Q 231 both assuming M = 0.55 Mp^; (3) e xpansion velocities neglect any 
momentum input other than radiation pressure on dust, based oniNenkova et ^ lll999li. 


mechanism pulsations can levit ate material from the star , 
enhancing the mass-loss rate (IWoodI 1 19791 : iBowe nl ll988ll . 
Strongly linked to this pulsation is the formation of dust 
in the denser ste l lar atmosph ere fe.g. lBladh &: Hofneill20i^ : 
[McDonald et ah] |2012| . |2014| ~). Radiation pressure on this 
dust drives the stellar wind in a cooler outflow, also of typ¬ 
ically ^1 0 km s~^, rising with stellar luminosity to ~20 km 
s“^ fe.g. ILoup et al.l[19931 : [Marshall et al.ll2004l 'l. Although 
dust production begins at luminosities b el ow th e RGB tip 
llSloan et al.ll201ol : [McDonald et ^I2011bl l3. l2012fl . sufficient 
radiation pressure to drive the wind i s only achieved once 
the star reaches ^1000-5000 Lo (e.g . Winters_eLal |[2000l: 
Riebel et al. | [2010[ : [ita fc Matsunagali201ll : [McDonald et alJ 


2014 ^ . 


However, the expansion velocities of metal-poor stel¬ 
lar winds are poorly known. The stars are smaller com¬ 
pared to metal-rich stars at the same luminosity, so pulsa - 
tions are thought to be weaker ([Kieldsen fc Bedding|[l99^ . 
The stellar surface resides deeper in the gravitational po¬ 
tential, so material must be levitated further before reach¬ 
ing temperatures where dust can condense. We also expect 
the dust:gas ratio will decrease, meaning that the associated 
mass-loss rate and velocity enhancements will also decrease. 
The transition between a chromospheric and pulsation- 
enhanced, dust-drive n wind should consequently happen at 
a higher lum inosity ([Winters et al.[ [20031 : [McDonald et al.) 
[2011al . [ 2 OI 2 II . Without measurement of these expansion ve¬ 
locities, however, these theories remain untested. The wind 
expansion velocity is an important parameter in deter¬ 
mining mass-loss rates from radiative transfer modelling 
of infrared excesses, which is a valuable method of de¬ 
termining the mass of both dust and gas ejected by gi¬ 


ant stars, and t he only method currently feasible outside 


our Galaxy (e.g. 

McDonald et al.[[2011al: [Bover et al.[|2011 

2 OI 2 I: Sloan et al 

2012[: Javadi et al.[2013[: Jones et ahlboiS 

Bover et al.l 2015 

)■ 


Once material leaves the star it is reprocessed, first by 
the interstellar radiation field (ISRF) and then by interstel¬ 
lar shocks. Other unknowns then become important in de¬ 
termining how quickly dust is destroyed, molecules are disso¬ 
ciated, and atoms ionised. Contributing factors include the 
star-formation history of the population (setting the UV flux 
from high-mass main-sequence, low-mass horizontal-branch 
and post-AGB stars), the amount of interstellar dust shield¬ 
ing, and the local stellar density. There is a growing real¬ 
isation that the enrivonment in which stars lie can have 

J irofound effects on the ir molecular and mineralogical yields 
Zhukovska et ^[2015h . 

Nowhere is this more strongly felt than globular clus¬ 
ters. The preparation of this work sparked the theory that 
ionisation should exert a powerful con trol over the interstel¬ 
lar en vironment of globular clusters ([McDonald fc Ziilstra[ 
[20153). Very strong ISRFs should be generated by post- 
AGB stars and the hottest white dwarfs, but the short-lived 
nature of these objects means the ISRF is constantly vary¬ 
ing. Not only should these sources be capable of dissociating 
and ionising the intra-cluster medium (ICM), but they also 
provide the ICM with enough thermal energy that it over¬ 
flows the globular cluster, preventing further star formation. 

We now turn our attention to the circumstellar envi¬ 
ronments of globular cluster stars. Many of the above issues 
can be addressed by observing (sub-)mm CO lines emanat¬ 
ing from the winds of stars in globular clusters. CO is one 
of the first molecules to form as the extended atmosphere of 
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the star is ejected, and it survives until its dissociation by 
the ISRF. This dissociation is mainly by photo-absorption 
in molecular lines, hence even in the absence of dust, self¬ 
shielding can be important in protecting CO in stellar winds. 
Observationally, the widths of (sub-)mm CO lines provide 
a direct measurement of the wind-expansion velocity, while 
their strength provides an independent estimate of mass-loss 
rate. Their strengths are also affected by the characteristic 
radius from the star at which CO is dissociated by the ISRF 
dMamon et aLlflQSSl l. Hence, by comparing the strength of 
CO lines from stars with known mass-loss rates to model 
predictions, we can quantify the strength of the ISRF to 
which they are subjected. 

In this paper, we report on observations with the Ata¬ 
cama Large Millimetre Array (ALMA) of the four brightest 
stars in 47 Tucanae. Section [5] describes these observations, 
and defines upper limits to the non-detections we find. Sec¬ 
tion [3] models the expected mass-loss rate from the star, the 
size of the CO envelope, and calculate the expected strength 
of the CO transitions. This is done under the assumption 
that the ISRF is strong dissociating the CO envelope. Re¬ 
sults are discussed in Section 2] Alternative scenarios which 
might remove the CO envelopes of our target stars are ex¬ 
plored in the Appendix. 


2 OBSERVATIONS 


2.1 Target selection 

Among the closest and most populous globular clusters is 
the intermediate-metallicity cluster 47 Tuc ( NGC 104; 4.5 
kpc, JFe/H1 = -0.72 dex, M = 1.5 x 10® Mm: lGnedin et ^ 
I 2 OO 2 I: iHarrid l20ldL This cluster was chosen partly be¬ 
cause its stars show a relatively small spread in chem- 


size llWorlev et al.l 2010l: 

Gratton et al.l 20131: Gordero et al.l 

2 OI 4 I: ICerniauskas et al. 

I 2 OI 4 I: iDobrovolskas et al.l l2014l: 

Lebzelter et al.l 2014: Thveesen et al.l 20141: Johnson et al.l 

2014|). and partly because it contains a well-studied set of 


luminous, dusty AGB stars. 

The cluster’s var iable stars have a unique, long and 
well-covered history dPickerind Il894l f . Infrared photome¬ 
try and spectra h ave been published at many wavelengths 
over many epochs dG IassfcJbasjilQTj; Origlia et al.lll99'79 : 


iRamdani fc JorissenI 200ll : l0riglia et al 2002). Mid-infrared 

spectroscopy of the brightest variables jLe bzelter et al.l 

I 2 OO 6 I : Ivan Loon et al.l l2006al : iMcDonald et al.l 2011al l show 
a mixture of silicate and alumina dust features, and a po- 
tential featureless con tribution from metallic iron dust (cf. 
iMcDonald et al.]|2010l l. This combination ha s led to sev¬ 
eral estimates of their dust-production ra t es dOriglia et al 
200^ BowretalJ 2010l:l0riglia et al.ll2010l: IMcDonald et al 


2011b Momanv et al.l[2012l L High-resolution optical spec¬ 


troscopy of moderately bright giants suggest relatively slow 
outflows from the stellar chr omospheres of ~10 km s“^ 
I MuDonjJd_j^TOnjmonl l2007ll . Near-infrared spectroscopy 
i Lebzelter et al.l 2014 1 shows that the brightest stars have 
abundances similar to other stars in the cluster. They have 
not experienced substantial third dredge-up, and perhaps 
have not experienced any third dredge-up episodes. 

We present sub-mm CO observations of the four most 


Table 3. Properties of the observed stellar winds derived from 
the upper limit to the CO(2—>1) line flux. 


ID 

CO(2-^l) 
(mjy kms“^) 

R.M.S. Noise 
(mJy kms“^) 

a 

VI 

<183 

85 

2.1 

V2 

<168 

114 

1.4 

V3 

<170 

99 

1.7 

V8 

<162 

70 

2.3 


Notes: Based on CO (2—>-1) line flux. The root-mean squared 
(r.m.s.) noise is the quadrature-summed average flux of a re¬ 
gion of identical width, offset by -1-100 km s“^ from the ‘de¬ 
tection’. The significance (cr) is the multiple of the r.m.s. noise 
by which the line ‘detection’ is above the mean. 


lu minous, most e volved stars in 47 Tuc. They are denoted 
by H awver HoggI (I 1973 II as VI, V2, V3 and V8. Previously 
published parameters of these stars from the above refer¬ 
ences are given in Tables [T] and [2] In particular, Table[2]gives 
the mass-loss rate (Mir) and wind velocity ( 1 ;®.)*°*) p r edicted 
for a purely radiation-driven wind. IMcDonald et al.) (l2flllah 
derived these on the basis of radiative transf er modelling of 
their infrared spectra using the dusty code (|Nenkova et al.l 
Il99!j l with a radiation-driven wind (density type = 3). In 
all four cases, the spectral energy distributions were best 
modelled with a mostly metallic iron wind, with some stars 
having contributions from amorphous and crystalline sili¬ 
cates, aluminium oxide and iron oxide. Table[2]also gives the 
mass-loss rate expected for a purely magneto-acoustically 
driven wind, as derived from semi-empirical scaling laws ap¬ 
plied to the cluster’s other stars (Mneimers, Mscos; see Sec¬ 
tion [32|. 


2.2 ALMA & APEX observations 

Observations of 47 Tuc VI, V2, V3 and V8 (see Figure [T|) 
were carried out using ALMA in receiver band 6 on the 
night of 2013 Nov 05 over a 1.05 h integration. Observations 
were centred on the ^^C^®0 J = 2 —>■ 1 transition at 230.538 
GHz, and were carried out at high precipitable water vapour 
(2.8-3.5 mm) but otherwise good conditions. The correlator 
output consisted of 3840 channels of 244.141 kHz each. The 
data were processed following the standard ALMA quality 
assurance calibration and imaging processes. Neptune was 
used for amplitude calibration and the quasars J1924-291 
and J0102-7546 for bandpass and phase calibration, respec¬ 
tively. 

During the ALMA observations the majority of anten¬ 
nae (26 of 29) were positioned within ~250m of the ar¬ 
ray centre; the remaining three telescopes were at larger 
distances (~530m, ~790m and ~1014m). At the imaging 
stage of data reduction the visibilities from these three an¬ 
tenna were downweightecj^. This was accomplished using the 
uvtaper parameter set to True a nd the outaper option set 
to 1 . 5 ” in CASA’s clean task llMcMullin et al.|l2007l ). In 
addition to this a restoring beam of 2" x 2" was also ap¬ 
plied at the GLEANing stage. These techniques were used 


^ This provides a trade-off between competing noise sources: 
noise is increased by the larger solid angle within the beam, but 
additional noise from the poorer phase calibration on longer base¬ 
lines is reduced. 
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Right Ascension 

Figure 1. Map of imaged regions (squares), with the cluster 
centre (green + symbol) and pulsars (small blue letters, from 
iFreire et alj |2003h shown. The large, green asterisk marks the 
only significant UV source: an optically bright post-AGB star. 


to generate images with a beam more representative of the 
antennae in the centre of the array. Imaging additional to 
that provided with the ALMA data release was conducted 
in CASA using the clean routine to create two data cubes 
for each source, one with a binning of ten spectral channels 
(channel width: 3.2 km s“^ bin~^) and one with two-channel 
binning (0.64 km s“^ bin“^). The pixel scale in both cubes 
is 0.15" pixel~^. 

No line or continuum source is visible in any of the maps 
to the depth of the observations. The standard deviation 
of the noise in the final images is ~14 mjy beam“^ per 
0.64 km s“^ channel pair for VI, V3 and V8, and ~16 mJy 
beam“^ per channel pair for V2. The shortest baseline (17 
m) corresponds to a maximum detectable object size of 15" 7, 
though the CO shells are expected to be unresolved (~l"; 
Section On-source spectra were extracted using a 3x3 
pixel box surrounding the central source. 

Figure [5] shows the two-channel-averaged spectra 
around the velocity of the cluster. The maximum two- 
channel flux for each source is 30, 33, 36 and 37 mJy beam“^ 
for VI, V2, V3 and V8, respectively. Each spectrum was ex¬ 
amined and found to have noise following a Gaussian dis¬ 
tribution, both on each source and at various points away 
from each source. On-source standard deviations in the spec¬ 
tra are 10.2, 11.9, 11.6 and 11.2 mJy beam“^ per 0.64 km 
s“^ bin for VI, V2, V3 and V8, respectively. 

Additional data were obtained from the Atacama 
Pathfinder Experiment (APEX) telescope. A four-hour in¬ 
tegration on 47 Tuc VI was made in the 345 GHz ^^C^®0 
J = 3 —>■ 2 transition (project identifier: 0-092.F-9327A). 
Standard reduction of the APEX data was performed at 
Onsala Space Observatory before receipt of the data. No 
line was detected, with a limiting peak flux of ~0.1 Jy over 
a typical 10 km s“^ bin from the 17'lZ beam. 



LSR velocity (km s'^) 

Figure 2. On-source ALMA spectra of the four observed sources 
(points), offset by an arbitrary flux level with z ero flux lines 
shown. The cluster mean velocity is -26.7 km s“^ llHarri ilili 
and the range of radial veloci ties over the pulsati on cycle are 
shown as thick horizontal lines llLebzelter et al .l200dl . Lines show 
the spectra, boxcar-smoothed by 11, 20, 15 and 7.4 km s“^ for 
VI, V2, V3 and V8, respectively. This smoothing corresponds to 
the upper mass-loss rate limits quoted in Table [3] 

2.3 Results & comparison to previous 
observations 

Neither the ALMA nor the APEX data showed any con¬ 
vincing indication of spectrally resolved CO emission (e.g. 
departure from Gaussian noise centred around zero flux). 

The only previous significant CO observations of the 
cluster are bv lOriglia et al.l ([l9973), who claim a weak CO 
(1—kO) detection around the cluster’s velocity with a signal- 
to-noise ratio of 3. We interpret their peak fluxes as 1.6 
Jy, assuming an aperture efficiency for the Swedish-ESO 
Sub-millimetre telescopcQ (SEST) of 27 Jy K“^. Origlia et 
al.’s observing position is not quoted, but assumed to be 
a « 00 24 22, 5 ~ —72 0 5 48, f ollowing their offset from 
iKrockenberger fc Grinding |T99^. This is 68" from VI, 
compared to a full-width half-maximum SEST beam width 
of 45". In neither dataset do we re cover any sou r ce com - 
parable to the GO (1-0) detection of lOriglia et al.l (|l997al L 
though the beam footprint of the SEST observations extends 
considerably beyond both the ALMA primary beam of 15"7 
and the APEX primary beam of 17"3. 

The CO emission is expected to be spatially unresolved 
by ALMA and the mass-loss rate is expected to be suffi¬ 
ciently low that the lines will be optically thin. An unre¬ 
solved, optically thin line will have a rectangular (‘boxcar’) 
spectral profile. To identify whether a low-contrast, spec¬ 
trally resolved feature is present in our data, we smoothed 
our data using a boxcar function (running average). The 

^ http://www.apex-telescope.org/sest/html/telescope- 
instruments / telescope / index, html. 
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boxcar was given different central velocities, from dlsr = 
—50 to 0 km s“^, and different widths, from 1.27 km s“^ to 
85 km s“^ (4-133 channels). 

Typical expansion velocities for dust-producing stars in 
the Milky Way are ~ 10-20 km s“^, and are typically ~10 km 
s~^ f o r stars producing little dust le.g. ISchoier &: OlofssonI 
I 2 OOII : iDanilovich et all l2015l h We presume the expansion 
velocity is Uexp < 20 km s“^, and that the centre of the 
emission line is Doppler shifted by no more than 5 km s“^ 
from the average radial velocity of the star. The limits to 
the velocity-integrated CO line strength are typically below 
<200 mjy km s“^, with a typical Icr uncertainty being ~100 
mjy km s“^ (see Table EJ. The boxcar-smoothed spectra 
producing these maximum fluxes are shown in Figure (2] 

The mean photospheric radial velocities and their 
pulsation-induced variation are listed in Table fffl. The high¬ 
est flux in the VI and V 8 spectra is within 10 km s“^ of 
th e mean photospheric v elocities of the stars as published 
by iLebzelter et al.l (l2005l f . These may indicate a tentative 
detection of the CO line, but we do not consider them reli¬ 
able as they are only at the 1.3 and 2.Icr confidence levels, 
respectively. 


3 THE EXPECTED CO LINE FLUX 

3.1 The CO flux under Galactic ISRF condtions 

Our observations were carried out under the naive assump¬ 
tion that the envelopes of our target stars behave as Calac- 
tic stars do. In this section, we show that we could expect a 
similar Galactic star to be observable at the distance of 47 
Tucanae. This provides a baseline to which we can compare 
our final model. 

I OlofssonI (I2OO8I I provides an estimate of the peak CO( 2 - 
1 ) line flux from an evolved star, namely: 



where M is the mass-loss rate in Mq yr“^, Ve is the expan¬ 
sion velocity in km s“^, fco is the ratio of CO to H 2 , and 
D is the distance in kpc (D = 4.5 kpc). From Table [21 we 
can modestly expect M « 2 x 10~^ — 2 x 10“® Mq yr~^, 
depending on the star and the method used. If we assume 
M > 3 X 10“^ Mq yr“^ (appropriate for VI), a canonical 
Ve = 10 km s“^ an d fco = 4 x 10 ~^ (ba sed on scaling a 
canonical 2 x 10“^ llRamstedt et al.ll2008lf by the metallic- 
ity of 47 Tuc, which is a fifth of solar metallicity), we obtain 
a peak flux of >14 mJy. Assuming a boxcar profile gives 
an integrated intensity of >280 mJy km s“^, which is still 
gre ater than o ur up per limit. 

lOlofssonI ([ 200^ 1 also predicts that the CO (3-2) flux 
should be around twice as strong. A less conservative value 
of M = 10“® Mq yr“^ produces 1200 mJy km s“^ for the 
CO(2-l) line, implying 2400 mJy km s“^ for the CO(3-2) 
line, or 1.2a in our APEX observation. 


® Note that the velocities published bv iLebzelter et ^ ll2005h are 
given as heliocentric velocities (T. Lebzelter, private communica¬ 
tion), although this is not explictly stated in the text. A correction 
of—8.7 km s“^ has been applied. 


Much depends on the exact values used in the equa¬ 
tion above, as well as the accuracy of the formula itself. 
Stars in globular clusters exist in environments unlike the 
well-studied stars in the Solar Neighbourhood. In particu¬ 
lar, it has recently be en shown that the ISRF is co nsiderably 
stronger and harder ([McDonald fc Ziilstrall2015^ ). This will 
have a substantial impact on the size of the CO envelope and 
the line flux emitted from it. In this section, we re-evaluate 
the expected CO line fluxes from these stars by creating a 
photochemical dissociation model for VI. We consider other 
factors peculiar to globular clusters in the Appendix, but 
find that they have little impact. 


3.2 Expected mass-loss rates 

3.2.1 Expected mass-loss rate from chromospheric physics 


In general, mass-loss rates of stars without dust-driven 
winds are comparatively well determined and have been pa- 
rameterised into scaling relations. Two common formalisms 
which are freq uently adopted are the emprical formalism of 
[Reimerd ([l97^) and its semi-emprical modification by SC05. 
In these formalisms, the mass-loss rate is scaled by an efh- 
ciency factor, 77 . In terms of observable units, they can be 
written as: 

Mr = 4 X (2) 


and: 


Msc = 4 X 10 ^®?7sc 



4300MTf‘/ ’ 


(3) 


where the luminosity, mass and temperature (L, M, T*) are 
all scaled to solar units and the mass-loss rate is in Mq 
yr“^. T hroughout a large sample of Milky Way globular 
clusters, [McDonald fc Ziilstra[ ([2015bh compared the mass 
lost between the main-sequence turn-off and the horizontal 
branch and compared it to theoretical isochrones at fixed 77 . 
This Axes rj around the point where mass loss is strongest, 
close to the RGB tip. From these, the following efficiency 
parameters were derived for the median stars in globular 
clusters: 


77 r = 0.477T 0.0701^62. and 

? 7 sc = 0.172 ± 0.024^H2^ (4) 


where each value is given with its respective statistical and 
systematic erroi0. From these, we derive the minimum mass- 
loss rates listed in Table [2] of Mueimers ~ 2 — 3 x 10~^ Mq 
yr“^ and Mscos ~ 7 — 11 x 10~^ Mq yr“^. These values 
are appropriate for ‘stable’ stars, without extra energy from 
pulsations or radiation pressure on dust. While the interplay 
of pulsations on transfer of magneto-acoustic energy to the 
chromosphere is poorly determined, we would expect these 
values to provide a lower limit to the mass-loss rate from 
pulsating, dusty AGB stars. 


iHevl et al.l 1 I 2 OI 5 I I model diffusion of AGB stars in 47 Tuc, ar¬ 
guing for a lower rj on the RGB (t/r ss 0.1), and a higher r] on the 
AGB (rjR ^ 0.7). This would increase the predicted mass-loss rate 
tor our targets. Private communication s with Heyl et al. have not 
identified the source of the discrepancy. [Lebzelter &: Wood] (I 2 OO 5 I I 
also argue for a slightly smaller ryR 0.33, on the basis of the 
stellar period-luminosity diagram. 
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3.2.2 Expected maximum rate from stellar evolutionary 
arguments 

The mass lost from an AGB star during the last part of 
its life must not exceed the envelope mass of the star at 
the beginning of that period. Low-luminosity AGB stars are 
not easily seperable from RGB stars. AGB stars can only be 
uniquely identified once they exceed the luminosity of the 
RGB tip. If we can determine both the remaining lifetime 
of an AGB star passing the luminosity of the RGB tip, and 
the envelope mass it has when it does so, the quotient of the 
two should yield the average mass-loss rate of AGB stars 
above the RGB-tip luminosity. 

The time spent above the RGB tip is relatively easy 
to compute. Our four stars lie above the RGB of 47 Tuc. 
For the purposes of these calculations, we also include the 
fifth-brightest star in the cluster. V4 is a probable AGB 
star at the luminosity of the RGB tip (see the Hertzsp rung- 
Russell diagram presented by [McDonald et al.ll2011ah . Like 
its brighter counterparts, it shows long-period pulsations 
and an infrared excess, consistent with circumstellar dust 
production. Assuming there are typically 5 ± stars 
above the cluster’s RGB tip, and taking an evolutionary 
rate of one star passing the RGB tip every 80 000 years 
^McDonald et all 1201 lal l. these five stars should represent 
the last 0.40 ± 0.18 Myr of AGB evolution in the cluster. 

The envelope mass at the luminosity of the RGB tip is 
more difficult to calculate. iGratton et al.l (|201fll l determine 
zero-age horizontal branch (ZAHB) star masses in 47 Tuc. 
Their ground-based measurements suggest the ZAHB star 
mass ranges from 0.629 M© to 0.666 Mq, with a median of 
0.648 Mq, but they also include masses derived from Hubble 
Space Telescope (HST) photometry, which range from 0.650 
Mq to 0.691 Mq, with a median of 0.674 Mq. We assume 
these stars will end their lives as a 0.53 Mq white dwarf, 
base d on measured ma s ses of white dwarfs i n other clus- 
ters (|Richer et al.l Il997l : iMoehler et al.l l2004l : iKalirai et al.l 
I 2 OO 9 I I. These masses imply that ~0.12 Mq or ~0.16 Mq is 
lost during the entire horizontal branch and AGB evolution, 
including that below the RGB-tip luminosity. Assuming the 
limiting case of losing 0.16 Mq over 221000 years, the av¬ 
erage mass-loss rate for stars above the RGB tip of 47 Tuc 
must be <7 x lO”’^ Mq yr“^. 

A more limiting case is reached if we can estimate the 
envelope mass at the luminosity of the RGB tip itself. Calcu¬ 
lating the mass lost between the zero-age horizontal branch 
up the AGB to the luminosity of the RGB tip is not triv¬ 
ial. Th e evolutionary track used by [McDonald fc Ziilstral 
|[2015b[ ) was designed to reproduce the RGB and HB of 47 
Tucanae. This track was generated with the MESA (Modules 
for E xperiments in Stellar As trophysics) stellar evolution 
code ([Paxton et akll^lll . [20l3 l. Mass loss was included at 
the rate of 77 R = 0.45, reproducing a ZAHB mass of 0.666 
Mq after 12.02 Gyr. The model loses a further 0.067 Mq be¬ 
tween the ZAHB and the point where the star reaches the 
same luminosity as the RGB tip (i.e. the same luminosity as 
V4). Taking the limiting case of the i/5T-derived mass of 
0.674 Mq, a total mass at the luminosity of the RGB tip of 
0.607 Mq is implied, giving an envelope mass of ~0.077 Mq 
and limiting the average mass-loss rate to <3.5 x 10“^ Mq 
yr“^. In contrast to the mass-loss rate in the previous para- 


Table 4. Modelled ISRF incident on each sta r at 1076 A , relative 
to the flux in the Solar Neighbourhood from [jur i <1974 162 000 
photons cm~^ s“^ A~^). The centile probabilities approximate 
the best estimate, 1-cr, 2-cr and full ranges of the modelled fluxes. 


Star 1076 A ISRF strength at given centiles 



0 

2.25 

15.85 

50 

84.15 

97.75 

100 

VI 

1.4 

1.8 

7.0 

49 

102 

170 

256 

V2 

0.55 

0.65 

3.6 

19 

37 

59 

87 

V3 

0.31 

0.47 

2.6 

9.0 

22 

124 

483 

V8 

4.3 

5.2 

14 

135 

315 

551 

858 


graph, this mass-loss rate requires that Reimers’ mass-loss 
formula is appropriate for early-AGB stars. 


3.2.3 Summary 

Reimers’ scaling lav0 provides a lower limit to the mass-loss 
rates for pulsating AGB stars, implying an average mass- 
loss rate for our four stars of >2.3 ± 1.1 x 10~ Mq yr“^. 
Stellar evolution models, pinned by the mass of HB stars 
f Section 13. 2. 211 suggest an upper limit of <3.5 x 10“^ Mq 
yr“^. We can therefore expect the mass-loss rate from our 
four observed stars to lie in the range 1.2-3.5 x 10“^ Mq 
yr"A 

This predicted mass-loss rate is estimated to produce 
a CO line flux within the noise limit of our observations. 
However, the average predicted mass-loss rate is much less 
than the average of 1.4 x 10“® Mq yr“^ derived from these 
stars’ infrared spectr41, although the infrared-based rates 
do come with a sizeable and poorly quantifiable uncertainty. 
These infrared-based mass-loss rates assumed slow (2-4 km 
s~^) winds, accelerated only by radiation pressure on dust; 
faster (>10 km s“^) winds would increase the calculated 
mass-loss rate. 


3.3 The interstellar radiation field (ISRF) 

The radiation environment within gl obular clusters differs 
from t hat in the Solar Neighbourhood ([McDonald fc Ziilstra[ 
[2015ah . The ionised hydrogen and lack of interstellar dust 
within globular clusters reduces the opacity for ionising pho¬ 
tons. Ionising sources are mainly post-AGB stars and white 
dwarfs, rather than the O and B stars of the Solar Neigh¬ 
bourhood. This leads to a much harsher ISRF than is typi¬ 
cally assumed. 

The ISRF of 47 Tuc was modelled in detail by 
[McDonald fc Ziilstra[ ([2015a[ L The ISRF can vary consid¬ 
erably, as the presence of very hot post-AGB stars in the 
cluster is stochastic. Many of the hottest post-AGB stars 
and white dwarfs in the cluster remain unobserved. We must 

® This process can be repeated with the SC05 scaling law, which 
predicts an average of >8.5 ± 4.3 X 10“^ Mg yr“^ for our stars. 
This is inconsistent with the <3.5 X 10“^ Mq yr“^ provided 
by the evolutionary analysis. This suggests the SC05 law is not 
appropriate for this particular situation. 

® Note that these assume a dust-to-gas ratio of 1:1076. The metal- 
licity of 47 Tuc limits the abundance of refractory elements such 
that t he true dust-to-gas ratio should not be closer to unity than 
1:778 ([McDonald et al1[2011ah . 
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Figure 3. The ISRF estimated for each star. The line denotes 
the best-estimate, time-averaged value, with lightening colours of 
grey representing the 68.3, 95.5 and 100th centil e ranges. The 
ISRF in the Solar Neighbourhood (as assumed bv iMamon et ahl 
Il988h is shown as the dotted green line. 

therefore rely on stellar evolution models to describe the typ¬ 
ical ISRF within 47 Tuc and its temporal range. 

To compute the ISRF impinging on each AGB star, 
we have created a three-dimensi onal model of the cluster . 
We begin with the catalogue from [McDonald et al.l (l2011lj) . 
which contains the two-dimensional (RA-Dec) position of 
almost every post-main-sequence source in the cluster up to 
the early post-AGB phase. We assume this is also represen¬ 
tative of the distribution of more-evolved, hotter post-AGB 
stars. For each catalogued source, we assign a physical depth 


within the cluster. Assuming an Earth-to-cluster distance of 
4500 pc, we randomly select a source nearby in right as¬ 
cension to each star and use the declination of the random 
source to assign a radial distance from the cluster centre for 
each target source. This creates a stellar distribution with 
the same radial density profile as the density profile in decli¬ 
nation, approximating a spherically symmetric distribution. 

We then take th e post-A GB stellar evolution model of 
[McDonald fc Ziilstral ([2015a[ l. following the emitted spec¬ 
trum as it evolves through its upper-AGB, post-AGB and 
hot-white-dwarf evolution. We follow the evolution of each 
of our model sources in time, randomly forcing one of our 
sources to undergo this post-AGB evolution with an average 
of onc e per 80 000 years (the stellar evolution rate accord¬ 
ing to [McDonaId et all[2011bl 'l. Summing the modelled flux 
from these sources at the locations of our four target stars, 
we obtain the integrated interstellar radiation field imping¬ 
ing upon them. 

This model was run 100 times, in order to simulate dif¬ 
ferent radial distributions of the radiation sources and ob¬ 
served targets. For each run, the evolution was followed for 
10 ® years to quantify the time variability of the radiation. 

Figure 0 shows the ISRF, averaged for each source over 
all runs and all times. For comparison, the figure also shows 
the Solar Neighbourhood ISRF used bv [Mamon et al.l ([l988[ . 
originally from [Jur^[l974[ i. The radiation field in the 914- 
1120 A region (normally associated wi th CO dissocia t ion) i s 
typically many times higher than the [Mamon et aP ([l988[ l 
value, and it is also considerably harsher. 

3.4 Modelling the outer CO shell radius 
3.4-1 A crude estimate 

CO in the circumstellar environment has some self-shielding, 
but most shi elding comes f r om ci rcumstellar dust and H 2 , 
as shown bv [Mamon et al] l[l988[ . their figure 2). At some 
depth the atmosphere becomes optically thick to incoming 
or outgoing radiation, forming the photosphere, or r = 1 
layer. In mass-losing stars, particularly at short wavelengths, 
the r = 1 layer is wavelength dependent, depending on the 
shielding within the wind. Dissociation can therefore take 
place in an unshielded environment, if the time to intercept 
a UV photon is large compared to the time taken to reach 
the r = 1 layer, or a shielded environment if it is smaller. 

If shielding is not important, the CO dissociation ra¬ 
dius will be limited simply by the average time it takes a 
CO molecule to encounter an interstellar UV photon, and 
the CO dissociation radius around metal-poor stars should 
be very similar to metal-rich stars. If the shielding is impor¬ 
tant, CO will be dissociated close to the r = 1 layer, and 
the dissociation radius around metal-poor stars should be 
smaller. However, it cannot become much smaller before H 2 
becomes the dominant shielding mechanism. In either case, 
we expect metallicity to only have a second-order effect on 
the CO dissociation radius, and that the primary effect will 
be the strength of the ISRF. For a Solar-Neighbourhood 
ISRF, we would expect the 47 Tuc stars have CO envelopes 
filled to ~0.013-0.023 pc (^Section lAlll . A stronger ISRF will 
dissociate CO proportionally more quickly, so the dissocia¬ 
tion radius should be inversely proportional to the ISRF 
strength. 
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Figure 4. Radial properties of the toy photo-chemistry model. 
This model does not include heating by the central star, hence the 
difference between the temperature profiles in the bottom panel. 


The primary contrib ution to CO dissoc iation is nor¬ 
mally the 1076 A band (|Mamon et al■lfl^88^ ■ though the 
harder radiation field may change this in our case. The po¬ 
tential range in ISRF strength is large (Table |4|) and is con¬ 
siderably greater for the stars closer to the cluster centre 
(VI, V8) than stars on the cluster’s outskirts (V2, V3). 
Taking the Icr ranges for the potential variation in ISRF, 
we compute a typical range of CO envelope outer radii of 
~4 X 10“® - 0.008 pc, or 8-1800 AU. This does not for 
changes in the (self-)shielding of CO or the harder ISRF, 
and implies CO may be dissociated very close to the stel¬ 
lar surface (<10 R*) at certain times. If the wind remains 
optically thin below these heights, there may be substantial 
CO dissociation within the dust-formation shell (2-10 R*) 
around these stars. More typically, however, we would ex¬ 
pect 90 per cent of the CO to be dissociated by ~100 AU 
(~120 R*) in VI, and ~30 AU (^40 R*) for V8. 


3.4.2 A simple model 

To get a simple estimate of the abundance of CO that could 
be expected in a high-UV-field, metal-poor environment we 


set up the following toy model, which w e base on 47 Tu c VI. 
We used the moving-mesh code AREPO llSpringelll201Cll l. Our 
version of the AREPO code has been modified to in c lude t ime- 
dependent che mistry as prese nted by Snfith_et_^ ll2014h ff or 
full details see IClove 3 l2007h and I Glover fc Clarkl ll2ol2 ')~). 
The CO abundance is calculated assuming that the CO for¬ 
mation is limited by an initial radiative association step, and 
that the CO destruction rate is primarily through photodis¬ 
sociation. The gas is shielded from the ambient radiation 
field by dust an d gas self-shi e lding^ which we calculate using 
the approach of IClark et al.l ll2012h . 

The physical model consists of an expanding shell of gas 
with a total mass of 9 x 10“^ Mq, following a M(r) oc 
power law, starting from the assumed stellar surface at R, = 
1.3 X 10^® cm = 0.87 AU. The entire envelope corresponds 
to the mass ejected over 3000 yiQ at a mass-loss rate of 3 
X 10“^ Mq yr y. The gas initially follows a T{r) oc 
power law (e.g. ICroeneweaenl l2012|j with the inner surface 
corresponding to the temperature of the stellar photosphere, 
taken to be 3517 K. 

In this oxygen-rich environment, the fractional abun¬ 
dance of CO is dete rmined by the carbon abundance: [C/Fe] 
= -0.13 ± 0.20 dex llRoedieer et al.ir2014h . The carbon abun¬ 
dance varies between [C/Fe] ~ -0.5 dex for a few CNO- 
processed stars to [C /Fe] = 0.0 dex fo r the majority of 
the stellar populationj[BrjIeY_et al.l l2004h . We use [Fe/H] = 
—0.72 dex, following lHarrisI (l2010h . and [C/H] = -0.72 dex to 
conservatively model the carbon-rich end o f the population. 
We as sume a dust-to-gas ratio of 1:1000 (iMcDonalcl et aP 
l2011ah . and [O/H] = -0.57 dex dRoediger et al.ll2014h . The 
envelope is initially fully molecular, with all its carbon in 
the form of CO. 

The CO abundance relative to protons was tracked just 
inside the inner boundary of the expanding shell, represent¬ 
ing an expanding region of gas in a dynamic envelope. The 
code assumes an external UV field with a spe ctrum equiva - 
lent to that seen in the Solar Neighbourhood llDrainelllll78l l. 
but with a flux that is 50x higher (i.e. 50 Habing unit^; see 
Section [3] and Table 0. 

This model includes several important simplifications, 
which generally lead to an overestimate of the CO abun¬ 
dance, and an overestimate of the CO dissociation radius. 
The most important of these assumptions are: 


• The adoption of an ISRF from the Solar Neighbour¬ 
hood, scaled upwards by a factor of 50 to the expected time- 
averaged flux of the ISRF in 47 Tuc at 1076 A. The UV field 
in globular clusters also contains harder UV photons, which 
will lead to more rapidly dissociated CO, and is likely to 
result in significant quantities of CO^, or even CO^"*". 

• A lack of radiative heating from the central star, and 
a lack of thermal heating by circumstellar dust (collisional 
heating is included). This means the temperature of the CO 
envelope is colder than we expect (Figure 0, protecting the 
CO in lower ro-vibrational states. Figure 2 of iMamon et al.l 


^ This value was chosen to be sufficiently larger than the typical 
envelope filling time of similar stars under solar-neighbourhood 
conditions (~1300-2200 yr; Section lAll l. The exact radius and 
time are not critical to the results. 

® A Habing unit is the integrated flux between 912 and 1110 A 
in the Solar Neighbourhood, i.e. 1.6 x 10“^ erg s“^ cm“^ 
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Table 5. Parameters used in the CO line flux models 


3.5 A SMMOL 


Parameter 

Symbol 

Unit 

Value 

Used in both the AREPO +photochemical model and SMMOL 

Gas temperature^^l 



1 

c 

Dust:gas ratio^^’^l 



ITOOO 

Initial CO/H 2 



5.33 x 10-® 

Used in the arepo +photochemical model 

Outflow velocity 


km s ^ 

10 

Wind density at R* 


cm“^ 

2.63 x 10® 

Wind density law 


cm“® 

(^)”" 

Mass-loss rate 

... Mq yr 1 

3 X 10-’^ 

Dust condensation 

-^cond 

cm 

1.23 X IQi® 

Dust temp, at -Rcond 


K 

[4] 

ISRF 


Habing 

50 


Used in SMMOL 


Stellar radius PI 

R* 

cm 

1.23 X IQi® 

Stellar spectrum 



[2] 

Stellar mass 

M, 

M0 

0.55 

Stellar luminosity 

L* 

L© 

4824 

stellar Teff 

T, 

K 

3623 

Dust condensation^^] 

-^cond 

cm 

2.46 X IQi® 

Dust temp, at -Rcond 


K 

850P1 

Dust composition 



[2] 


Notes: [1] The photo-chemical model does not include the cen¬ 
tral star so the gas cools radiatively. Heat exchange with dust 
is implemented, but dust and gas heating by the central star 
are not, so the ejecta cools radiatively. Smmol resets the tem- 
perature structu r e to t his assumed power law. [2] Following 
[McDonald et al] [3] Based on scaling a solar value of 

1:200 by [Fe/H] = -0.72 dex. The ratio is limited by condensable 
metals to <1:778 llMcDonald et al]|2011a(j . [4] As parameterised 
by kinetic temperature law. 


Smmol is a spherical-geometry, non-LTE radiative trans¬ 
fer code which solves level populations with the ac- 
cel erated lambda iter a tion (ALI) scheme, described 
by Rvbicki &; Hummerl (|l99lh and IScharmer &; CarlssonI 
(|l985[l . It relies on a discrete spatial grid and ray tracing 
for radiation transfer and is optimised for high-optical-depth 
regimes. 

Sm mol adopts dust properties from [McDonald et ah] 
( 2Qllal), which uses the Dusty radiative transfer code 
I Nenkova et al ]|I^ i,o model the spectral energy distri¬ 
bution of VI. This model uses a gr ain mixture of 88 per 
cent metallic iron dOrdal et al.lll988l l and 12 per cent sili¬ 
cates toraine fc Ledl 19841 '). surrounding a star modelled by a 
BT-SETTL model atmosphere spectrum (I Allard et al.ll2003l ~l. 
The gas temperature structure was taken from our photo¬ 
chemical modelling results. An accelerating wind was used, 
with a velocity at the inner region of 8 kms“^, and a ter¬ 
minal velocity of 11 kms“^, to mimic radiative acceleration 
of dust. This lowers the optical depth of the proHle but, as 
our stars have low mass-loss rates, the resulting line flux dif¬ 
fers negligibly from a constant 10 km s“^ wind. A turbulent 
velocity of 2kms“^ was adopted. Table [5] lists a full set of 
physical parameters. 

Smmol modelled the observation as observed by a single 
250-metre telescope with a perfect Gaussian beam. The CO 
envelope (Cf.'05) remains a point source at this resolution 
(l”07), such that the conversion from antenna temperature 
to flux units (Janskys) represents the flux from the entire 
source as observed by ALMA. 


(|l988h shows that halving the CO excitation temperature 
has a small but noticeable effect on CO shell size. 

Figure |4] shows the abundance of CO relative to pro¬ 
tons as a function of dynamic radius. At the stellar surface, 
the carbon is entirely bound in molecules, with a CO:proton 
ratio of 2.67 x 10“®. The abundance drops rapidly with dis¬ 
tance due to photodissociation by the ambient UV field. By 
a radius of 2.8 x 10^® cm (0.00092 pc, or 218 R*), 90 per 
cent of the CO has been dissociat ed. This repres e nts a shell 
~22x smaller than predicted by iMamon et al.l (Il988f) for 
a similar Galactic star. Bearing in mind the simplihcations 
above, we can expect the shell size to be smaller still, sug¬ 
gesting a close to inverse relation between ISRF strength 
and CO shell size. 

3.5 Predicting the CO line flnxes 

Millimetre/sub- mm CO lines at these m ass-loss rates are op¬ 
tically thin (e.g. iRamstedt et al.ll2008l L Given the r~^ den¬ 
sity law of our wind, we can expect the CO line flux to 
scale with the shell size, such that the CO line flux should 
be ~22x lower than predicted by IRamstedt et al.l (l2008l l 
for a star losing mass at 3 x 10“^ Mq yr“^. Other fac¬ 
tors are important, notably the temperature structure of 
the atmosphere, which controls the rotational level popu¬ 
lations of CO. We have therefore predicted line strengths 
from our envelope using the SMMOL radiative transfer code 
llRawlings fc Yat^l2001^ . Table [5] lists the parameters used 
in the model. 


3.5.2 Results 

Table [5] lists the resulting line fluxes from SMMOL. Smmol 
clearly predicts line strengths below the upper limits of our 
observations, showing that photo-dissociation explains well 
our non-detection of the CO emission from these stars. 

We stress that these fluxes are only order-of-magnitude 
estimates, as many wind parameters are not accurately 
known. The line flux is strongly affected by the tempera¬ 
ture structure of the inner stellar wind. We have adopted a 
simple power law for our model, but the physical tempera¬ 
ture will be signihcantly affected by the dust properties. The 
dust properties of our target stars are not well known, partic¬ 
ularly the grain mineralogy, grain size and dust-production 
rate, though they are likely very different from AGB stars 
in the Solar Neighbourhood llMcDonald et al] 1201 fal l. This 
uncertainty therefore translates into a considerable uncer¬ 
tainty on the CO line strengths. When combined with the 
time-variability of the ISRF (Section 13.311 . which can ex¬ 
ceed an order of magnitude, it becomes clear that we are 
limited to stating that the line fluxes will scale very approx¬ 
imately with ISRF strength, and that we would typically 
expect them to be below the sensitivity limit of our obser¬ 
vations. 


4 DISCUSSION & CONCLUSIONS 

In this paper, we have presented new ALMA observations of 
the four brightest stars in the globular cluster 47 Tucanae, 
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Table 6. Predicted line fluxes from the photo-chemical model of 
VI. Note these are conservatively bright estimates. 


Line 

Observational 

limit 

(mjy kms“^) 

Predicted line fluxes 
from SMMOL 
(mjy km s“^) 

J = 4-:>3 


16 

J = 3-i>2 

<2000W 

6.1 

J = 2^1 

<183 ± 8512] 

2.6 

J = l-:>0 


0.1 

^APEX, ^ALMA. See Section 12.31 


examining the CO J = 2 —>■ 1 transitions. No source was 
detected. Having modelled the interaction between the ISRF 
and the circumstellar enviroment, we expect that the CO 
line fluxes from the giant stars in 47 Tuc are around two 
orders of magnitude below our observation limit. 

For any given globular cluster giant, if the wind is opti¬ 
cally thin, then the radius of CO dissociation should exhibit 
an approximately linear relation with metallicity and an ap¬ 
proximately inverse relation with ISRF strength (Section[3]). 
The radiation field should be highly stochastic, driven by 
the temporally variable population of hot white dwarfs. The 
population of young white dwarfs is not well determined 
in most clusters, therefore it is not immediately possible to 
identify which clusters will be the best candidates for ob¬ 
servation. However, the ISRF can be expected to vary by 
several orders of magnitude within a cluster, both tempo¬ 
rally and spatially (Tabled)). Targets far from the cluster 
centre may experience radiation fields that are weaker than 
those in the Solar Neighbourhood (though they will still con¬ 
tain harder radiation). Giant stars at large distances from 
their host clusters’ centres may therefore provide the best 
metal-poor stars around which to detect stellar winds and 
measure their expansion velocities. Table|T]indicates that V2 
and V3 may occasionally experience ISRFs that are weaker 
than the Solar Neighbourhood. While this seems not to be 
the case at present, we can hope that a similar situation can 
be found in other Milky Way globular clusters. 

Since the higher stellar density forces the CO envelopes 
to be smaller than their Galactic counterparts, the CO(3“ 
2 ) line (or higher rotational states) may be considerably 
brighter than CO (2-1), as these are better populated in the 
warmer regions close to the star (Table [S]). We therefore 
suggest observations to detect CO line strengths and widths 
from globular cluster stars focus on: 

(i) stars at large radii from their host clusters and 

(ii) high CO rotational states emitted closer to the star. 


We also advocate further observations of the cir¬ 
cumstellar dust. Unexpectedly large am ounts of dust are 
seen around many globular c luster stars dSloan et al.ll201^ : 
[McDonald et al.ll2010l. 2011al'l. Their spectra have strong sil¬ 
icate features dU ebzelte r et al.l [20061 : Ivan Loon et al.ll20()6al : 
[McDonald et aLl2011al '). but also a n underlying continuum , 
possibly explained by metallic iron ([McDonald et aHl2010ll . 
This contrasts w ith the notable ab s ence of interstellar dust 
within clusters (|BpYer_et_aI. 200^ Matemag a_et_alj_ 200^ 


MtdDonald et al.l 20091 : Hover et al.[ 200ii : Barmbv et all 

2009l 'l. and suggests dust is being destroyed while still in 


the circumstellar environment. The highest-energy photons 
(~40 eV, ~300 A) are likely to penetrate into the dust¬ 


forming regions. While they are of insufficient number to 
dissociate a large fraction of CO in this region, they may 
still be important in the photo-chemical formation of dust. 

Finally, we note that this situation is unlikely to be 
limited to globular clusters. The dissociation of CO around 
AGB stars by nearby post-AGB stars and white dwarfs is 
largely a function of stellar density (although stellar age, 
elemental composition and the density of the surrounding 
ISM are also important). This source of CO dissociation 
may also become important in old open clusters, nuclear 
star clusters, and gas-poor galaxies. The implications for 
the formation, chemistry, evolution and survival of dust in 
such environments are significant. 
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APPENDIX A: OTHER PROCESSES 
AFFECTING THE STELLAR ENVELOPE 

In the main text, we have shown that photo-dissociation of 
CO by interstellar UV radiation is a significant mechanism 
which reduces the CO envelope size. In this Appendix, we 
consider whether alternative processes acting on globular 
cluster stars may also affect the size of the CO envelope. 


Al Envelope size, filling time and resolution 

In determining the effect of other processes, the relevant size- 
scales and time-scales must be compared to the extent and 
filling time expected for the CO envelope. In this section, 
we determine what those would be for an AGB star with 
the properties of those in 47 Tuc but with a Galactic ISRF 
incident upon it. We also examine the possibility that the 
CO might be resolved by the ALMA interferometer. 

The limiting size of a CO envelope in a Galactic en- 
vironment is usuall y set by its dissociation by the ISRF. 
[Mamon et al.[ ([l988[ . their figure 3) predict that 90 per cent 
of the CO in the envelope of a Galactic AGB star with a 
mass-loss rate of a few x 10“^ Mq yr“^ will be dissociated 
by a radius of ~0.013-0.023 pc. This equates to ~0.6”-l.l” 
at the distance of 47 Tuc. Assuming a typical wind veloc- 
ity for less-evolved stars of 10 km s“^ ([Winters et al.[[2003[ : 
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lGroenewegen|[20l3 ) . a typical CO envelope would therefore 
be filled within 1300 to 2200 years. If CO is more rapidly 
dissociated by interstellar UV, we can expect the radius to 
be smaller and the filling time to be correspondingly shorter. 

This would preclude a much larger envelope, which 
could be resolved by the ALMA beam. Since the winds are 
optically thin, the dissociation radius is mainly set by the 
time taken for a CO molecule to absorb a UV photon. A 
faster wind may therefore lead to a larger envelope, but the 
velocity and envelope size will not be directly proportional. 
Larger envelopes will be less self-shielding, as the optical 
depth in the dissociating lines will decrease. It is unlikely 
that an envelope could reach the velocity needed to extend 
beyond the 2” ALMA beam (>20 km s“^), and virtually im¬ 
possible for it to reach the >60 km s“^ needed to be diluted 
by a sufficient factor that it becomes unobservable in the 
central beam. This would still be detectable by the larger 
15.7” X 15.7” box, so we consider it very unlikely that the 
large-scale structure of our envelopes are resolved by the 
ALMA interferometer. 


A2 Episodic or variable mass loss 

If mass loss is episodic or variable, it is feasible that we 
could see strong mid-IR emission from recent dust produc¬ 
tion near the star, while the CO envelope further from the 
star traces a period of lower mass-loss rate. Civen the chro¬ 
mospheric and CO mass-loss rates are discrepant by a factor 
of ~10 in each case, any variability or episodic nature must 
be correspondingly large, i.e. a mass-loss rate variation of 
10 X or, alternatively, mass loss occurring only 10 per cent 
of the time. While the chances of finding any four stars going 
through this state is high, the chances of finding four partic¬ 
ular stars (in this case the most luminous) simultaneously 
going through an episode of enhanced mass loss is therefore 
^lO'^:! against. 

This probability only holds if the enhancement in dust 
production is uncorrelated with stellar luminosity. Such a 
correlation occu rs, for example, duri ng the bright phase of 
a thermal pulse. ICirardi et al.l (l20ld 'l model the thermally- 
pulsating (TP) phase of a globular-cluster-like ACB star to 
be ~1.2-1.8 Myr in length. If a star leaves the ACB every 
80 000 years jMcDonald et all l2011bh , this lifetime equates 
to 15-23 stars on the TP-ACB, four of which we observe, the 
remainder of which will be at luminosities at or below the 
RCB tip. The bright phase of a thermal pulse is r elatively 
short compared to the pulse cycle: ~0.3 per cent llHerwid 
I 2 OO 5 I L The chances of a particular star being in the bright 
phase of a thermal pulse are therefore ~0.3 per cent, and 
the chances for four such stars are ~1 in 12 billion. 


A3 Self-destruction by chromospheric heating 
and/or stellar UV 

Although C O bands are present in the near-IR spectra of 
these stars (iLebzelter et al.ll2005l L it is possible that they 
could produce enough UV photons to destroy their own CO 
envelopes. Metal-poor stars are warmer than their solar- 
metallicity counterparts due to reduced metal-line opacity 
in their atmospheres. Escaping stellar UV could lead to the 
destruction of CO near the stellar surface. A BT-Settl 


model atmosphere jAlIard et al.ll2003l at 4000 Lq and 3300 
K emits 3 x 10^® photons per second below the 1120-A CO 
photo-dissociation limit. By comparison, a mass-loss rate of 
3.4 X 10“^ Mq yr“^ and efficient condensation (i.e. a CO:H 2 
mass ratio of 1:500) yields 3 x 10"^^ CO molecules produced 
per second. Many more CO molecules are produced than 
there are UV photons to dissociate them. 

While the BT-Settl atmospheres do not model chro¬ 
mospheric or shock heating in the outer atmosphere, the 
heating cannot compensate for the factor of 10^® between 
the modelled number of ionising photons and number re¬ 
quired to dissociate all the CO the star produces. Assum¬ 
ing a typical wavelength of 1105 A, 3 x 10^^ photons pro¬ 
vide an energy of 8 x 10^® W, or 10“^“^ W m~^ as observed 
from Earth. Escaping photons would therefore be easily de¬ 
tectable in the GALEX FUV filter for all realistic spectral 
energy distributions. None of the four sources is present in 
catalogues from the satellit^, which is complete over 47 Tuc 
to a far-UV AB magnitude of «20.0 («1 x 10~^® W m“^), 
approximately two orders of magnitude more sensitive than 
would be necessary to detect such a UV field. 

Reduced metal-line and dust cooling in pulsation-driven 
shocks presumably mean that the post-shocked gas further 
out in the atmosphere also cools more slowly, leading to more 
effective dissociation of CO further from the stellar surface. 
However, observations show the expected column depths of 
CO absorption lines in the H-, K^- and ( archival) L-band 
VLT/CRIRES data llLebzelter et al.ll2005l L indicating that 
CO must survive to at least a few stellar radii. The observed 
CO absorption and lack of observed UV emission from these 
stars means that self-destruction of CO by the stars that 
produce it appears unlikely. 


A4 Stellar encounters 


Stellar encounters are common in globular clusters, and are 
thought to be responsible for a lack of pla nets and binary 
stars , and a prevalence of s tellar e xotica (lGiUiland__et_alJ 
2 OOOI: [Pavies fc Sieurdss^ I 2 OOI : IWeldrake et al.l l2005l : 


Knigge et al.l 20081 : Snurzem et al, |2009| L While stellar en¬ 
counter rates have previously been computed for 47 Tuc, 
published valu es are not specifically a pplicable to giant stars’ 
envelopes (e.g. iDavies fc Benj[i995l L 

The encounter timescale (renc) between an ACB star 
and any other star in the cluster can be computed as: 


7 X 10^° yr 


10® pc-® 


Rp) Mp, 


10 kms ® R enc <M*>’ 


(Al) 


where n is the stellar number density, Venc is the typical rel¬ 
ative velocity at the start of the encounter, Renc is the maxi¬ 
mum distance at which an encounter is deemed to occur and 
< Mt > is the typical stellar mass (e.g. iDavies fc SigurdssonI 
I 2 OOIII . We conservatively assume an envelope radius of 0.03 
pc ISection lAlll . hence conservatively adopt Renc = 0.06 pc. 
We assume Venc = 22 km s”®, or twice the radi al componen t 
of the central velocity dispersion of the cluster (iHarri 32 OIOIL 
Civen an average stellar density within th e half-light radius 
of 1720 Mq pc”® and dN/dm oc m~^'^ jMarks fc Kroup^ 


http://galex.stsci.edu/ 
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l2010h . bounded by a planetary mass (0.08 Mq) and the max¬ 
imum sub-giant mass (~0.8 Mq), we obtain < M* > > 0.2 
Mq and n < 8600 pc“®. 

We conservatively derive Tenc > 340 000 years, though 
the possible range of values allow for Tenc > 1 Myr. This is 
much longer than the 2300-year envelope filling timescale. 
While a fraction of stars may undergo envelope-disrupting 
stellar encounters at some point in their dust-producing 
phase, this is extremely unlikely to affect all four observed 
stars simultaneously. 


A5 Disruption in the cluster potential 

A circumstellar envelope will experience tidal disruption or 
‘spaghettification’ as it interacts with the cluster poten¬ 
tial, which becomes significant once it experiences signifi¬ 
cant differential acceleration over the filling timescale (which 
is considerably shorter than the Myr orbital period). 
We can assume a 0.06-pc wide envelope in a radial orbit in 
alPlummeil (Il91ll i potential of characteristic radius 7.8 pc 
dLane et al.ll201(]f ). and calculate the change in distance be¬ 
tween opposite ends of the envelope. We find a maximum 
differential acceleration of 2 x 10“^^ m s“^, corresponding 
to a tidal elongation of one part in 10“*. Tidal deformation 
of the envelope is therefore negligible. 


A6 Sweeping by hot Halo gas 

There is a gener al paucity of interstellar media in globu- 
lar clusters ie.g. ISrnith et al1ll99d : [van Loon et al.lf2006bl . 
I 2 OO 9 II. with interstellar media detected in only M15 and 47 
Tuc (iFreire et al. gOOlJ: l^ans et al]l2003l : ISover et al]l2006l : 
Ivan Loon et al.ll2006br) . This strongly suggests that globu¬ 
lar clusters are efficiently clear ed of their ICM b y some pro¬ 
cess, on timescales of <1 Myr dSover et al.ll200^ 'l. We model 
that this is because it is ionised by hot white dwarfs and 
post-AGB stars, before being later cleared on slightl y longer 
timescales (a few Myr; (iMcDonald fc Ziilstrall^lSal ll. How¬ 
ever, this model does not include ablation by the passage of 
the cluster through the hot gas of the su rrounding Galactic 
Halo , which is another leading contender dFaulkner &: SmfM 

[TiiS). 

The stand-off distance of a bow shock (Ro), i.e. the 
minimum distance from the bow shock to an isolated, mass¬ 
losing star, depends on four parameters: the mass-loss rate of 
the star and the expansion velocity of its wind (M, Uexp), the 
interstellar medium density (pism) and the relative velocity 
of the star with r espect to the Halo gas (t^. Ro is given by 
dCox et al.l d2012ll . and references thereirFm: 


Ro 


MVe: 


47rpisMu; 


(A2) 


Depending on the rotational coupling between the 
Halo and Milky Way, u* lies between 57 km s“^ (co- 
rotating Halo) and 190 km s~^ (non-rotating Halo) 
dKrockenberger fc Grindlavlll995l i . The internal velocity dis¬ 
persion within 47 Tuc is negligible in comparison (11 km 


Note that this momentum balance assumes cool Halo gas. Bal¬ 
ancing ram pressure against gas pressure gives values of Rq about 
three times larger. 


s- MHarris| [2oIo|). No good determinations exist for pism; 
estimates at 3.1 kpc from the Galactic Plane are typi¬ 
cally nn ~ 0.01 but vary by an o rder of magni¬ 

tude d Albert et al.lll99^ . iTavlor fc CordesI dl993f l estimate 
an electron density of tie ~ 0.007 cm“® at this location (we 
can expect rie « nn in th e ionised Halo gas) , and t he more- 
recent spherical model bv iMiller fc BregmanI d2013ll predicts 
a lower value of rie « 0.0002-0.002 cm“^. 

Assuming M « 10“’’ — 10“® Mq yr“’^, Uexp = 4-15 
km s“’^ and nn ~ 10“^ — 10“’^ cm“®, Ro — 0.02-38 pc, 
always greater than the expected CO shell radius. Addi¬ 
tionally, these stars ar e not isolated, but are surr ounded 
by a hot, ionised ICM dMcDonald fc Ziilstra|[^15a^ . which 
should provide additional shielding out to Ro ~ 0.3-98 pc. 
We conclude that ram-pressure stripping by Halo gas should 
not be effective at removing circumstellar CO envelopes. 
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